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We consider analytically the dynamics of an arbitrary number and configuration of vortices in an 
annular Bose-Einstein condensate obtaining expressions for the free energy and vortex precession 
rates to logarithmic accuracy. We also obtain lower bounds for the lifetime of a single vortex in 
the annulus. Our results enable a closed-form analytic treatment of vortex-vortex interactions in 
the annulus that is exact in the incompressible limit. The incompressible hydrodynamics that is 
developed here paves the way for more general analytical treatments of vortex dynamics in non- 
simply connected geometries. 
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I. INTRODUCTION 

The experimental observation of quantised persistent 
currents in toroidally confined Bose-Einstein condensates 
(BECs) [IHS] is closely related to the appearance of vor¬ 
tices with quantised circulation HIS]. The dynamics of 
quantised vortices and dark solitons [SHE] in BECs is a 
key ingredient in the understanding of the properties of 
superconductors [9] and quantum liquids [10], and have 
been intensely studied by the BEC community [ni[in. 

The persistent current state is protected by an energy 
gap, preventing the creation of excitations at low ener¬ 
gies. Above the Landau critical velocity, which equals 
to the speed of sound in a homogeneous condensate, the 
superfluid state becomes unstable against the formation 
of excitations such as vortices m and dark solitons m 
in the bulk. The behaviour of these collective non-linear 
excitations can be distinctively different in non-simply 
connected geometries nsnn] as a result of the ‘holed’ 
topology. For example, in the ring geometry, the outer 
and inner surfaces play important separate roles in the 
vortex nucleation above the critical rotation [mUEI. The 
ring geometries are within experimental reach [IMI],for 
example, in a genuinely non-trivial ring topology where 
the central region is occupied by a tapered optical fibre 
and therefore essentially at infinite potential jUj. The 
ring traps can be highly flexible with a tunable radius 
and transverse frequency [2d] , making them ideal for the 
study of multiply-connected BECs. 

In a ring trap sustaining a persistent current, a tun¬ 
able weak link across one side of the torus has been 
used to shed vortices [IIMIIIS]. In another experiment, 
an annular condensate was stirred with a narrow blue- 
detuned optical beam resulting in vortices appearing in 
the bulk [^ . The weak link alters the local critical veloc¬ 
ity affecting the current in the ring, and forms a closed- 
loop atom circuit with applications in atomtronics |26| . 
In annular geometry, many-vortex states such as a circu¬ 
lar array of vortex-antivortex pairs can also be produced 


via the snake instability of a ring dark soliton p7I129j . 
whose subsequent vortex dynamics has been studied in a 
circular container [301 ED, and in a harmonic plus quartic 
trap EH EE] . It was found that when the BEC is rotated 
sufficiently rapidly, an annular geometry with a vortex 
lattice that evolves into a ring of vortices can emerge [33] . 

In a harmonic trap, an off-axis vortex precesses in the 
same sense as the vortex circulation, following a circu¬ 
lar trajectory [4]. This effect can be explained in terms 
of a Magnus force [53IEB] that is pointing radially in¬ 
wards. The Magnus effect is a transverse force that al¬ 
ways appears when a vortex moves with respect to a su¬ 
perfluid m, which has been linked to the concept of a 
geometric Berry phase [38] . In an infinite homogeneous 
neutral superfluid at zero temperature (in the absence of 
any quasi-particle scattering [39]), the transverse Magnus 
force relies on the superfluidity and single-valuedness of 
the condensate order parameter [35] ■ Vortex precession 
is a key observable feature of vortex dynamics, and here 
we seek to provide closed-form analytical results for vor¬ 
tex dynamics in the ring geometry, including precession 
in a ring trap. 

Another approach for studying vortex dynamics is 
formed by the method of images [iniiini mi, where 
we assume that the spatial vortex separation is much 
larger than the healing length. Replacing the conden¬ 
sate boundaries with image vortices provides a direct 
method for describing the solution to the hydrodynamic 
boundary value problem in the incompressible Thomas- 
Fermi limit. The method of images has been studied 
numerically for an annulus [HlllllllE]; here, we develop 
closed-form analytical expressions for key observable fea¬ 
tures describing an arbitrary number and configuration 
of vortices in a ring BEC that solve the boundary value 
problem exactly. We obtain expressions for the free en¬ 
ergy and vortex precession rates to logarithmic accuracy. 
Our energy agrees with that obtained by Fetter [44] in 
studying the critical rotation for exciting vortices in an 
annulus. On the other hand, our expression for the vor- 
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tex precession rate, an exact solution in the Thomas- 
Fermi limit, is qualitatively similar to previously numeri¬ 
cally and analytically approximately obtained precession 
rates [12]. As an application, these results are used to 
obtain a lower bound for the lifetime of a single vortex 
in a ring-shaped condensate. Furthermore, we show how 
the Uniqueness Theorem implies that the Magnus force 
approach gives equivalent results with the method of im¬ 
ages for a non-dissipative system in a flat-bottom box 
trap. 

The aim of this work is to elucidate observable aspects 
of vortex dynamics in an annulus in terms of closed- 
form results valid for an incompressible two-dimensional 
superfluid. Starting from the incompressible hydrody¬ 
namics presented here, it is possible to generalise the re¬ 
sults to spin-orbit coupled and compressible superfluids 
in non-simply connected geometries. The presentation is 
organised as follows: In Section |II A[ we derive an ex¬ 
pression for the free energy of an arbitrary number and 


configuration of vortices in an annulus. In Section IIB 


we apply the results to derive a closed-form expression 
for the precession rate of a vortex in an annulus. Unlike 
in a cylindrical container, the vortex can precess in ei¬ 
ther clockwise or counter-clockwise sense depending on 
its radial displacement. In Section |ll C we derive lower 
bounds for the lifetime of a single vortex in a ring trap. 
We defer technical details to the Appendices. 
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FIG. 1. (a) The free energy of a vortex in a uniform annular 
condensate, Eq. ®, with no = 1, ^ = 0.01, a = 0.5, and 
b = 4.0 (solid lin^ Also the Thomas-Fermi result is shown 
(dotted line) with rtrap = (a -|- b)/2. (b) The correspond¬ 

ing precession rates around the annulus, Eqs. 0 and ( [To| |. 
There is a quasi-stable radius at which w = 0 corresponding 
to the stationary point of the energy. The dashed-dotted lines 
correspond to the case as a —>■ 0 of a uniform condensate. 


energy due to the presence of the vortex) at position p 
(radial displacement of p) in a uniform annular BEG of 
density no occupying the area U is then given by (see 
Appendix 


II. VORTEX IN AN ANNULUS 
A. Free energy 

We wish to obtain the free energy E of having Ny 
vortices at {pj} {j = 1 , 2 ,... Ny) in the annulus fl = {a < 
r < b;0 < (p < 27r}. Throughout this work we scale our 
units so that h = 2m = 1, which means we measure time, 
length, and energy in terms of Gosc = 
and hujx respectively, where uj^ is the angular frequency 
of the trap in the a;-direction, and m is the mass of the 
atoms in the cloud. First focussing on a single vortex at 
p G U, we will work in terms of a stream function Xp(r) 
such that the superfluid velocity v is given by 

v(r) = —zx Vxp(r). (1) 

m 

It follows by definition of the vorticity V x v of irrota- 
tional flow that the stream function satisfies the Poisson 
equation. To take into account the physical no-flow re¬ 
quirement nov • n = 0 and the geometry of the system, 
where uq is the superfluid density and h the unit outward 
normal, we require Dirichlet boundary conditions: 

\ X{p,} (a) =X{p,-}(&) = «. 

We solve the hydrodynamic boundary value problem 
using the method of images (Appendix]^. The free en¬ 
ergy E = 2no |Vxpp of a vortex (i.e. the additional 
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= In 


47rno 

where ^ is the healing length, and 

>”<Uto(s) 
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^ 2n6^"(piP2)"(&^" — 


cos (n^ia), 
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where the bar denotes complex conjugation. All the 
geometry-related information is in the function g. We 
have defined 771 = pie“^b 772 = P 2 e“^^, {pi,P 2 } S 
[a, &];{ 0 i,^ 2 } G [0,27r[, and (j)i 2 = (j)i - (j) 2 . Equa¬ 
tion is given to logarithmic accuracy, i.e. we assume 
In ( 6 /^ ^ 1 , and that terms of order 0{^^) can be ne¬ 
glected. In the limit as a —)■ 0, Eq. ([^ simplifies to the 
well-known result [SIIISI 


lim E = dTTTio 

a—>-0 

for a uniform circular condensate of radius b with an off- 
axis vortex at r = p. 

For a uniform annular condensate with Ny vortices at 
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FIG. 2. Trajectories of a single vortex, (a) Free energy con¬ 
tours of two equal vortices such that one vortex is held fixed 
at r = 7 and 0 = 27r/3 (*) with a = 0.5 and 6 = 10 [see 
Eq. 1^]. (b) Same as (a) but with opposite vortices, (c) 12 
alternating vortices with 11 held hxed in a necklace of radius 
9 (a = 8 and b = 10). (d) Same as (c) but all the vortices are 
equal. The red thick curves show the condensate boundary. 


{pj} of charges {kj} respectively, we obtain 


E 

iirno 


Q{Pk < Pi) 

0(pfc > pi) 






= In ( - j + '^g{pk,Pk) + '^KkKi { 
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\Pk-Pi\ 

b 
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( 6 ) 


where 0 is the Heaviside step function, and /pki = 4>k — (t>i 
is the angle between the fcth and Zth vortex. Equation ([^ 
agrees with the result obtained by Fetter [44] in terms of 
the Jacobi theta functions; here, the physical interpre¬ 
tation is transparent and follows directly from the sepa¬ 
rately given contributions of the annular geometry (func¬ 
tions g) and the non-geometry related vortex terms. 

In the Thomas-Fermi approximation, valid for large 
strongly-repulsive condensates {b/^ » 1), we assume 
that the spatial variation of the density is so small that 
the kinetic energy comes only from the velocity, i.e. gra¬ 
dient of the phase. In particular, in this work we assume 
that Vno ~ 0, and the density at the vortex can be re¬ 
placed by the vortex-free state, as corrections are of order 
and the vortex does not alter the condensate 
density significantly. Denoting a disk of radius ^ centered 
at the vortex at p hy the free energy depends only 
on surface terms dB^ if we ignore the volume term aris¬ 
ing from Vno 7^ 0- Therefore, we need only multiply the 


free energy by the Thomas-Fermi density nTF(/5) at the 
vortex: 


E 


,TF 



f P- ’’trap 

V b 



(7) 


We have plotted the free energy of a single vortex [see 
Eq. 0] in Fig. 0a). The assumptions underlying an 
incompressible superfluid break down when the vortex is 
within a few healing lengths from the borders. 


B. Vortex precession in an annulns 


As the Gross-Pitaevskii equation is Hamiltonian in 
structure, at T = 0 in absence of dissipation, a vortex will 
follow equipotential contours of the free energy, which 
in general depend on density inhomogeneities. We have 
plotted contours of the V^-vortex free energy in Fig. 
They are the trajectories of a single (positive) vortex if 
the other vortices are held fixed. In general, the iV„ vor¬ 
tices will follow simultaneous trajectories that leave E 
time-invariant. 

As the free energy of a single vortex depends only on 
the radial coordinate, the cylindrical symmetry means 
the vortex will precess around the trap centre, and the 
rate can be calculated using the Magnus force. 


nm{K X p) = \^pE, 


( 8 ) 


where n is the number density and k = {hK/m)z is the 
circulation vector of a vortex of charge k G Z. Using 
Eqs. 0 and 0, in a uniform condensate in the annulus 
D, an off-axis vortex at r = p will precess at the rate 



while in an annular Thomas-Fermi condensate nxF = 
no 1 — () , the precession rate is 


wtf = uj -\- 
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( 10 ) 


We show plots of Eqs. 0 and (10) in Fig. 0b). It 
can be seen in Fig. 0b) that the annular boundary can 
formally stabilise a vortex (i.e. w = 0), which is not the 
case in cylindrical geometry. The precession frequency 
changes sign at this quasi-stable radius. Similar results 
have been obtained in Ref. |42] using numerical integra¬ 
tion of the Gross-Pitaevskii equation and the method 
of images, and an approximative analytical solution of 
the Poisson equation, where the annulus is conformally 
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mapped onto a straight strip. In comparison, our solu¬ 
tion for the stream function is an exact Green’s function 
of the Dirichlet problem in the annulus. 

Instead of the Magnus force approach of Eq. Q, the 
precession around the ring could have been obtained by 
numerically solving equations of motion arising from the 
method of images, updating the most current image con¬ 
figuration at every time step [H]. At first sight it may 
seem strange that the Magnus force and the method of 
images give equivalent results, as any equivalence might 
be expected to be ruled out because the method of im¬ 
ages concerns the boundaries of the condensate, while 
the Magnus force is a property of the vortex itself [SS] 
present even in an infinite homogeneous system. How¬ 
ever, the boundaries affect the free energy of a vortex 
(the functions g), and when there is no dissipation, the 
vortex must follow contours of the free energy. There¬ 
fore, when the vortex moves, the transverse Magnus force 
must become exactly the force arising from the gradient 
of the free energy, and thus be affected by the boundaries. 
The method of images, on the other hand, works because 
by the Uniqueness Theorem any solution to the Poisson 
equation with given boundary conditions [see Eq. ([^] is 
unique. For an incompressible superfluid in a flat-bottom 
trap, therefore, the free energy contours can thus be ob¬ 
tained either with the method of images, or by consider¬ 
ing the Magnus force. 


C. Phonon radiation by a processing vortex 

In Ref. [46], the power P radiated by a vortex exe¬ 
cuting circular motion in an infinite uniform system was 
obtained to be 


4c2 ’ 


( 11 ) 


where = —h^2Trno/Tn, and the speed of sound c = 
yJgnQjm. Setting P = dE/At = {dE/dp)dp/dt and us¬ 
ing Eqs. @ and §, we may obtain a lower bound for 
the vortex lifetimes T^a,b} in the annulus: 


QE dp 

4c^ 

QE dp 


7rg2 

Jr^±wi dp w^p^ 


( 12 ) 


where is the radius corresponding to w = 0. Numer¬ 
ically evaluating the integral (Fig. shows that while 
the vortex lifetimes remain small for very narrow conden¬ 
sates, the vortex quickly becomes long-lived as the thick¬ 
ness of the annulus is increased. Decreasing the healing 
length also makes the vortex lifetime longer. 


III. CONCLUSIONS 



FIG. 3. Vortex lifetime, (a) Lower bound for the lifetime of 
a vortex with g = ±1 drifting towards the inner edge r = a. 
Here b = 10.0. (b) Same as (a) but for a vortex drifting 

towards the outer edge r — b. Here a = 0.5. In both figures 
5 = 0.01, and we have scaled c = 


rates, and vortex lifetimes. The results form an exact 
solution to the incompressible hydrodynamic boundary 
value problem, and provide closed-form expressions for 
key observable features of vortex dynamics in an annu¬ 
lus. The incompressible hydrodynamics presented here 
forms the starting point for more generalised discussions 
of vortex dynamics in spin-orbit coupled and compress¬ 
ible superfluids |47j in non-simply connected geometries, 
which will be discussed later. 
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Turku Centre for Quantum Physics for hospitality. 


Appendix A: Green’s fnnction in an annulus 

In this Appendix, we derive the Green’s function in 
an annulus ft = {a < r < b;0 < ip < 2 tt} with Dirichlet 
boundary conditions using the method of images. By def¬ 
inition, the Green’s function satisfies the Poisson equa¬ 
tion 


V^Gpir) = 27TS{r- p), (Al) 

where J is a two-dimensional delta function, and {r, p} € 
ft. In general, the Green’s function can be decomposed 
into a radially symmetric singular and a regular part, 
which we label by Ep{r) and Hp{r) respectively: 

Gp{r)=Fpir) + Hp{r). (A2) 

The singular part Ep{r) is the fundamental solution cor¬ 
responding to system response at a general point r to a 
point source at some arbitrary point p, and can be writ¬ 
ten down as 


We have studied vortex dynamics in an annular BEG, 
obtaining expressions for the free energy, precession 


Ppi^) = -^ln(|r-p|^) 


(A3) 
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The regular part Hp{r) is harmonic in n satisfying 
Laplace’s equation there. The choice of this analytic har¬ 
monic function is asserted by the boundary conditions. 
Installing N images at {pj}, all outside of fl, we can 
write 


i=i 


where the sign depends on the winding number of the 
particular image vortex. 

For the annulus, we need two sets of images that both 
have an infinite number of image vortices [221 ESI EH]- 
Importantly, the limit of infinitely many images must 
be taken with care at the origin when constructing the 
Green’s function below. The image sets are constructed 
by considering the outer and inner surfaces of the an¬ 
nulus as separate circular boundaries, whose effect on 
each other must be compensated by a limit of having 
two infinite series of images. Assuming a vortex with 
circulation n at p, the outer surface induces an image 
at Pi = {b‘^/p)p of circulation —k. While temporarily 
fixing the boundary condition at the outer surface, the 
inner surface must be compensated with another image 
at P 2 — {a^/b'^)p of circulation k, as well as an image 
vortex at the origin of circulation —k [49] . which then 
affect the outer surface again, and so on. A similar chain 
of images is needed starting from the inner surface. 

Considering a stream function Xp(r) such that the su¬ 
perfluid velocity v is given by 


v(r) 


n 

m 


z X Vxp(r), 


(A5) 


it follows by definition of the vorticity V x v of irrota- 
tional flow that the stream function satisfies the Poisson 
equation. To take into account the physical no-flow re¬ 
quirement nv • n = 0 and the geometry of the system, 
where n is the superfluid density and h the unit outward 
normal, we require Dirichlet boundary conditions. In the 
annulus fl with Ny vortices at {Pj} {j = 1,2,... Ny), we 
have: 


ditions, we obtain 




= - In 


00 - 


- |2 
^-Pj\ 




(9 


-2iVln(-l - 

.a 


r-p| / 

In (2) In (6)-f In (^) In (r) 




i=i ^ 
+ ln 


/O' 

\ 2 J 52 , 

2 " 


1 —P 


V 5/ 

P 



— (a y-)- 6) 




-f (a -o- &) / . 


(A7) 


Equation (A7) for Xp(r) can be regrouped to read 


Xp(r) = 0(r > p) 
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cos [n{(p -(/>)] I , 
(A8) 


where 0 = re“^, 77 = pe"^ is the position of the vortex, 
and 0 is the Heaviside step function. We show a few 
example plots of Eq. (A8) in Fig. 

In the limit as a —>■ 0 of a circular trap, the stream 
function in Eq. (A81 is given by 


lim Xp(r) = In 

a—>-0 


152 - rpe‘(‘^-'^)| 

b\j.Qli‘p-'t’) — p| 


(A9) 


which can be decomposed into the singular and regular 
parts In (5/|re4‘^“‘^) — p|) and In (|l — /b‘^\), re¬ 

spectively. 


Appendix B: Evaluation of the free energy 


The free energy of a vortex at position p in a uniform 

{ j^ annular BEG of density no occupying the area H = {a < 

^ X{pj }(!■) = 27r^^.^i 5(r — Pj ), r < 5; 0 < (p < 27 r} is: 

X{pa(“)=X{pa(^) = 0- 1 f 

E = -no / |7;|2 = 00 . (Bl) 

^ Jn 

The diverging part arises from the vortex core, and does 
Inserting the images as per the Dirichlet boundary con- not contain physically interesting information. Let us 
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FIG. 4. Annular BEG 
The stream function x I 


with a = 0.5, and b = 4.0. (a) 

.e. constant-velocity flow lines, see 
Eq. (A8|] of a single vortex (*) at p = 2.0. (b) Two (positive) 
vortices at p = 2.0 separated by an angle of 27r/9. (c) Same 
as (b) but with two opposite vortices. The red circles show 
the condensate boundaries. 


where dn denotes the outward unit normal derivative on 
dB^. Using Eq. (A8), the individual terms are readily 
evaluated to yield 

Inf 


E = Anno { In ( - | -b 


In 


(p2-a2)(&2_p2) 


— a^) 


(B3) 


+ 2 E 


°° „2nfu2n 
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_ p2n'^(^p2n _ ^2n^ 

2n6^"p2"(6^” — a^") 


therefore regulate the kinetic energy by excising i?|, a 
small disk of radius ^ centered at the vortex, where 
^ is the healing length. Using Eqs. ( |A5[ ) and (A2), 
V • (nyVy) = (Vn) • (yV^) + n|Vxp -b nxV^y, = 
27r(5(r — p), and the Divergence Theorem noting that the 
singularity is not inside the integration domain, the finite 


This result is given to logarithmic accuracy, i.e. we as¬ 
sume In {h/^) 1, and that terms of order can be 

neglected. 

Having multiple vortices in the condensate only adds 
non-singular interaction terms of the form 

which can be evaluated using similar methods as 
above. For a uniform annular condensate with Ny vor¬ 
tices at p = {pj} of charges {kj} respectively, we obtain 
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(B4) 


where (j)ki = 4>k — 4>i is the angle between the /cth and Zth 
vortex. 

When a —)■ 0, Eq. (B4) simplifies to 
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